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A field-theoretical description of the behavior of homogeneous, elastically isotropic, compressible 
systems characterized by two order parameters at the bicritical and tetracritical points is presented. 
For three-dimensional Ising-like systems, a similar description is performed in the two-loop approx- 
imation in three dimensions. The renormalization group equations are analyzed, and fixed points 
corresponding to different types of multicritical behavior are determined. It is shown that the effect 
of elastic strains causes a change from a bicritical behavior to a tetracritical one and leads to the 
appearance of a wide variety of multicritical points. 
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Earlier, it was shown [1] that, because of striction ef- 
fects, elastic strains lead to the appearance of multicrit- 
ical points that are absent in the phase diagrams of cor- 
responding incompressible substances. 

The subject of this paper is the study of the influ- 
ence of striction effects on systems whose phase diagrams 
already contain multicritical points of the bicritical or 
tetracritical type. In the first multicritical point 

corresponds to the intersection of two lines of secondorder 
phase transitions and one line of first-order phase transi- 
tions, and in the second case, it corresponds to the inter- 
section of four lines of second-order phase transitions. In 
the immediate vicinity of a multicritical point, the sys- 
tem exhibits a specific critical behavior characterized by 
the competition between different types of ordering. As 
a result, at a bicritical point, one critical parameter is 
displaced by another. A tetracritical point allows the ex- 
istence of a mixed phase with the coexistence of different 
types of ordering. Such systems [2] can be described by 
introducing two order parameters that are transformed 
according to two irreducible representations. 

In structural phase transitions that occur in the ab- 
sence of the piezoelectric effect in the paraphase, elastic 
strains act as a secondary order parameter whose fluc- 
tuations are not critical in most cases [3]. Since, in the 
critical region, the main contribution to the striction ef- 
fects comes from the dependence of the exchange integral 
on the distance, only elastically isotropic systems are con- 
sidered in this paper. 

The model Hamiltonian of the system has the form 



H = I dPx 
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Here, <&(x) and ty(x) are the fluctuating order param- 
eters; Moi and uq2 are positive constants; t\ ~ \T — 
T c i\/T c i and T2 ~ \T— T c2 \/T c2 , where T c \ andT C 2 are the 
phase transition temperatures for the first- and second- 

3 

order parameters, respectively; y(x) = ^ u ott (i), where 

a=l 

u a fj is the strain tensor; g\ and g 2 are the quadratic stric- 
tion parameters; (3 is a constant characterizing the elas- 
tic properties of the crystal; and D is the space dimen- 
sion. In this Hamiltonian, integration with respect to 
the components that depend on the nonfluctuating vari- 
ables, which do not interact with the order parameters, 
has already been already performed. 

Changing to the Fourier transforms of the variables in 
Eq. (1), one obtains the Hamiltonian of the system in 
the form 

H Q = 1 J d D q{r 1 + q 2 )<P q <S>- q 
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Here, the components y$ describing uniform strains are 
separated. According to [1], such a separation is neces- 
sary, because the nonuniform strains y q are responsible 
for the acoustic phonon exchange and lead to long-range 
interactions, which are absent for uniform strains. 

Let us determine the effective Hamiltonian that de- 
pends on only the strongly fluctuating order parameters 
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$ and * of the system as follows: 

exp{-ff [*,*]} = B J exp{-# [$, ( 3 ) 

If the experiment is performed at constant volume, the 
quantity is a constant, and the integration in Eq. (3) 
only goes over the nonuniform strains, while the uniform 
strains do not contribute to the effective Hamiltonian. 
In an experiment at constant pressure, the term Pil. is 
added to the Hamiltonian, with the volume being repre- 
sented in terms of the strain tensor components in the 
form 



phase diagram in the fluctuation region are determined 
by the CallanSymanzik renormalization group equation 
for the vertex parts of the irreducible Greens functions. 
To calculate the f3 and 7 functions as functions involved 
in the CallanSymanzik equation for rcnormalized inter- 
action vertices m, 112, U3, gi, 52, <7i°\ 52°^ or complex ver- 
tices z\, Z2, w\, W2, Vi, V2, V3, which are more conve- 
nient for the determination of the multicritical behavior, 
a standard method based on the Feynman diagram tech- 
nique and on the renormalization procedure was used [6] . 
As a result, the following expressions were obtained for 
the functions in the two-loop approximation: 



ft = O [l + Uaa + u a a Uf3f3 + 0(u 3 )} (4) 
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and the integration in Eq. (3) also performed over 
the uniform strains. According to [4], the inclusion of 
quadratic terms in Eq. (4) may be important at high 
pressures and for crystals with strong striction effects. 
As a result, one obtains 



H= \ I d D q(r 1 +q 2 )^^ q 



q(T 2 +q 2 )^ q ^- 



i JdD 

+ lf / dD 1i(^ q l^q2)(^ q 3^- q l- q 2- q 3) 
P J d D ft($ 9 l$ 92 )(* 93 *- 9 l- (Z 2- 9 3) 

d D qi (<t> ql <t>_ ql )(<t> q2 <S>_ q2 ) 

d D q t (<f ql *_ ql )(* q2 *_ q2 ) 



(5) 



2^03 

4! 



Pvl 



-Vi 



n + 



m 



A 1 Z 

6~ Vl + 6 W3 



41n+ 190 , 23m 2 2m , 

Vi ViVn K, 

243 1 243 3 27 3 



Pv 1 



-v 2 



m + 8 2 n 



-7. v 2 + 

6 2 6 3 



41m + 190 
243 ' 



23n 



2n 



"243^- 27**' 
2 2 

Pv3 = ~V\ + -U3 + 



(6) 



(n + 2) 



W1W3 



m + 2 



5(n + m) + 72 3 23(n + 2) 



6 



-v 2 v 3 



486 



486 



23(m + 2) 



v 2 v 3 



n + 2 



-V1V3 



V1V3 
m + 2 



-V2V3, 



486 " 2 "° 9 

n + 2 ■ - 3 , o_„ -.2 171 



Pzi = -z\ H — V1Z1 + 2nz\ + 2mz 1 z 2 + —v 3 z 2 

23(n + 2) , 7m 2 2m 2 
243— ^ ~ 243^ " ¥ V2 ' 

m + 2 „ 3 „ 2 71 

Pz2 = -22 H 5 — v 2 z 2 + 2mzg + 2nz 1 z 2 + -v 3 zi 



23(m + 2) 



7n 



2n 



243 243^- 27^' 



(3 wl = -wi + 



n + 2 



V1W1 + \nz\w\ — 2mw\ 



+ (^ 2 ~ ™) Jd D q,(%^-^^- q 2) +Amz lZ 2W2 - 2m Wl w\ + ^v 3 w 2 



v 01 = u 01 - 12Zi, v a2 =u 02 ~12z 2 , 
VQ3 = u 03 - I2z 1 z 2 , Z 1 = -j=, 

.92 g° 



23(n + 2) 2 7m 2 2m 2 

vtwi viwi V0W2, 

243 1 243 3 27 3 ' 



•w2 



-W2 



m + 2 



V2W2 + Amz 2 w 2 — 2nw 2 



2 n 

+Anz\z 2 w\ — 2nw 1 W2 + —V3W1 — 

O 



w 2 



23(m + 2) 



In 



2n 



This Hamiltonian leads to a wide variety of multicriti- 
cal points. As for incompressible systems, both tetracrit- 
ical (v 3 + \2{z x z 2 - wiw 2 )) 2 < (vi + 12(zf - w\)){v2 + 
12(z 2 — w 2 )) and bicritical (u 3 + \2{z\z 2 — wiw 2 )) 2 > 
(v\ + 12(z 2 — w\))(v2 + 12(z| — w%)) behaviors are pos- 
sible. In addition, the striction effects may give rise to 
multicritical points of higher orders. 

In the framework of the field-thcorctical approach [5] , 
the asymptotic critical behavior and the structure of the 
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It is well known that the perturbative series expan- 
sions are asymptotic and the vertices of the interactions 
of the order parameter fluctuations in the fluctuation re- 
gion are sufficiently large for Eqs. (6) to be directly ap- 
plied. Therefore, to extract the necessary physical in- 
formation from the expressions derived above, the Pade- 
Borel method generalized to the multiparameter case was 
used. The corresponding direct and inverse Borel trans- 
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formations have the form 

f(vl,v2,v3,zl,z2,wl,w2) = 
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For the analytic continuation of the Borel transform 
of a function, a series in an auxiliary variable is intro- 
duced: 

F(vi, v 2 , v 3 ,z-i, z 2 , w 1 ,w 2 , 0) = (9) 

oo 
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and the [L/M] Pade approximation is applied to this 
series at the point 6 = 1. This approach was proposed 
and tested in [7] for describing the critical behavior of 
systems characterized by several vertices of interaction 
of the order-parameter fluctuations. The property of the 
system retaining its symmetry when using the Pade ap- 
proximants in variable 9 is essential for the description 
of multivertex models. 

In the two-loop approximation, the functions were cal- 
culated using the [2/1] approximant. The character of 
critical behavior is determined by the existence of a sta- 
ble fixed point satisfying the set of equations 

p i {v* 1 ,v* 2 ,v* 3 ,z* 1 ,z* 2 ,w* u w* 2 ) = Q {i = !,...,!). (10) 

The requirement that the fixed point be stable is re- 
duced to the condition that the eigenvalues bi of the ma- 
trix 

Bi ' j = ^ 

Vi,Vj = v 1 ,v 2 ,v 3 ,z 1 ,z 2 ,w 1 ,w 2 ) 



lie in the right complex half-plane. 

The resulting set of resummcd (3 functions contains a 
wide variety of fixed points lying in the physical region 
of the vertex values with vi Vi > 0. 

A complete analysis of the fixed points, each of them 
corresponding to the critical behavior of a single order 
parameter, was presented in our recent publication [8]. 
Now, I consider the combined critical behavior of two 
order parameters. 

The analysis of the values and stability of the fixed 
points offers a number of conclusions. The bicritical fixed 
point of an incompressible system (v\ = 0.934982, v 2 = 
0.934982, v 3 = 0.934982, z x = 0, z 2 = 0,wi = 0,w 2 = 0) 
is unstable under the effect of uniform strains (&i = 
0.090, b 2 = 0.523, b 3 = 0.667, 6 4 = 0.521, b 5 = 0.002, b 6 = 
0.521, b 7 = 0.002). The striction effects lead to the 
stabilization of the tetracritical fixed point of a com- 
pressible system (n = 0.934982, v 2 = 0.934982, v 3 = 
0.934982, zi = 0, z 2 = 0, wi = 0, w 2 = 0, h = 0.090, b 2 = 
0.523, 6 3 = 0.667, 64 = 2.144, b 5 = 0.267, b 6 = 5.223, b 7 = 
0.882). 

The stability of other multicritical points cannot be 
investigated in terms of the described model, because the 
calculations lead to a degenerate set of equations. The 
degeneracy is eliminated by considering the Hamiltonian 
with allowance for the terms of higher orders in both 
the strain tensor components and the fluctuating order 
parameters. 

Thus, the striction-caused interaction of the fluctuat- 
ing order parameters with elastic strains leads to the 
transition from the bicritical behavior to the tetracrit- 
ical one and also to the appearance of new multicritical 
points with their own types of critical behavior in the 
phase diagram of the substance. 
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